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ABSTRACT 


A  method  has  been  developed  to  predict  three-dimensional 
hypersonic  laminar  or  turbulent  shock-layer  flows  for  perfect 
gas  or  equilibrium  air.  A  two-layer  eddy-viscosity  model  is 
used  for  the  turbulent  regime.  The  thermodynamic  and  trans¬ 
port  properties  for  air  are  obtained  by  interpolation  within 
a  two-dimensional  table  or  from  curve-fit  data.  Comparisons 
are  made  for  air  in  chemical  equilibrium  and  perfect  gas  for 
a  seven-degree  half-angle  spherically  blunted  cone  at  various 
flight  altitudes  with  a  cold,  moderately  cool  and  an  adiabatic 
wall  for  angles  of  attack  up  to  twenty  degrees.  Wall  heat 
transfer,  wall  pressure,  force  and  moment  coefficients  and  ex¬ 
ecution  times  are  compared  for  some  sample  cases.  This  method 
can  be  used  to  predict  viscous  flow  in  chemical  equilibrium 
over  axisymmetric  reentry  vehicles  at  angles  of  attack  up  to 
25  degrees. 


INTRODUCTION 


y 


There  is  renewed  interest  in  the  problem  of  computing  the  hypersonic  flow 
past  a  blunt  body  under  flight  conditions.  Due  to  the  current  interest  in  the 
improved  accuracy  of  ballistic  and  lifting  reentry  vehicles,  a  perfect  gas 
model  can  no  longer  accurately  predict  the  thermodynamic  and  transport  proper¬ 
ties  required  for  Inviscid  and  viscous  flowfleld  analyses.  Earlier  investigators 
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used  a  table  look-up  procedure*  or  curve-fit  data^  to  model  the  thermodynamic 
and  transport  properties  of  air  in  chemical  equilibrium  for  viscous  boundary- 
layer  flows. 

Reentry  vehicles  operate  through  a  wide  range  of  flow  conditions.  The 
complex  flowfield  is  bounded  by  the  body  and  the  bow  shock.  Crossflow  sepa¬ 
ration  may  be  present,  and  the  viscous  effects  may  predominate  over  the  entire 
flowfield.  The  full  Navier-Stokes  equations  are  elliptic  in  all  three  space 
directions,  and  a  numerical  solution  is  difficult  and  requires  large  com¬ 
puting  times  and  storage.  For  moderately  high  values  of  Reynolds  number, 
there  is  no  need  to  solve  the  full  Navier-Stokes  equations.^  The  classical 
approach  of  dividing  the  flowfield  into  a  viscous  boundary-layer  region  and 
an  outer  inviscid  region  becomes  inaccurate  for  reentry  flowfields  due  to  the 
presence  of  viscous  effects  throughout  the  entire  flowfield.  The  parabolized 
Navier-Stokes  (PNS)  approach*  uses  a  parabolic  approximation  in  the  streamwise 
direction.  This  involves  large  matrix  solutions,  and  the  computing  times  are 
still  quite  large.  The  viscous  shock-layer  (VSL)  approach**  developed  by 
Murray  and  Lewis  for  three-dimensional  flows  is  parabolic  in  both  the  stream- 
wise  and  crossflow  directions.  Since  the  crossflow  momentum  equation  is  para¬ 
bolic,  the  crossflow  separated  region  on  the  leeward  side  cannot  be  treated. 

The  solution  for  the  windward  region  up  to  crossflow  separation  is  accurate, 
and  the  computing  times  are  relatively  small.  In  the  viscous  shock-layer  solu¬ 
tion,  the  entire  flowfield  from  the  body  to  the  shock  is  treated  with  a  uni¬ 
form  set  of  equations.  The  problems  associated  with  the  displacement  thick¬ 
ness  interaction  and  edge  conditions  are  eliminated,  with  vorticity  interaction 
present  in  the  inviscid  +  boundary-layer  approach.  As  a  result,  the  viscous 
shock-layer  method  treats  all  higher-order ■ boundary-layer  effects  (displacement, 
vorticity  interaction,  longitudinal  and  transverse  curvature,  including  proper 
matching  conditions)  in  a  straight-torward  and  consistent  manner,  making  the 
viscous  shock-layer  approach  especially  attractive  for  design  studies. 

Recently,  a  numerical  method  was  developed  to  predict  laminar,  transitional 
and/or  turbulent  hypersonic  flows  for  a  perfect  gas^  over  a  blunt  body  at  angle 
of  attack.  In  that  approach  a  two-layer  eddy-viscosity  model  proposed  by 
Cebeci®  and  the  transition  model  developed  by  Dhawan  and  Narasimha^  were  used. 
This  method  has  been  extended  to  include  the  effect  of  air  in  chemical  equi-- 
librium.  Results  from  a  table  look-up  procedure  and  curve-fit  methods  have 
been  compared  with  those  from  a  perfect  gas  analysis.  The  table  look-up  pro¬ 
cedure,  described  later,  is  known  to  model  the  thermodynamic  properties  of  air 
quite  accurately,  while  the  transport  properties  are  curve-fits  of  Hansen®  data. 

Numerical  solutions  are  presented  for  a  seven-degree  half-angle  spherically 
blunted  cone  at  a  flight  velocity  of  20,000  ft/sec  at  altitudes  of  40,000, 

80,000  and  160,000  ft  with  a  cold  (540R) ,  moderately  cool  (3600R)  and  an 
adiabatic  wall  for  angles  of  attack  up  to  twenty  degrees.  Wall  heat-transfer, 
wall  pressure,  force  and  moment  coefficients,  center  of  pressure  locations  and 
execution  times  are  compared  for  perfect  gas  and  air  in  chemical  equilibrium 
for  some  sample  cases. 


ANALYSIS 

The  basic,  three-dimensional,  viscous  shock-layer  equations  are  derived 
from  the  steady  Navier-Stokes  equations  in  a  surface-oriented  coordinate  system 
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(s,n,i^').  The  governing  equations  for  turbulent  flow  are  developed  using 
methods  analogous  to  those  presented  in  References  9  and  10.  The  normal  ve¬ 
locity  V  and  normal  coordinate  n  are  assumed  to  be  of  order  e  ,  and  second- 
order  terms  are  retained  in  the  s-momentum,  (); -momentum  and  energy  equations. 

The  nondimensional  turbulent -conservative  form  of  the  equations  in  a  body- 
oriented  coordinate  system  can  be  found  in  Reference  5. 

EQUATION  OF  STATE 

P  =  p(p»h) 

For  a  perfect  gas,  the  above  equation  has  the  analytical  form: 

p  =  yp/[(y-i)t] 

For  a  gas  in  chemical  equilibrium,  the  functional  relation  may  be  given 
by  a  table  or  an  approximating  analytical  expression  (curve-fit)  as  discussed 
in  the  section  on  thermodynamic  and  transport  properties. 

The  same  equations  are  valid  for  laminar  flow  if  the  turbulent  eddy- 
viscosity  e'*’  is  set  to  zero. 

For  perfect  gas  flows,  since  the  Prandtl  number  is  not  a  function  of  n, 
it  can  be  taken  out  of  the  derivative  in  the  energy  equation.  For  flows  in 
chemical  equilibrium,  the  Prandtl  number  is  a  function  of  pressure  and 
enthalpy  and  thus  of  n;  therefore,  it  must  be  retained  within  the  derivative 
term. 

BOUNDARY  CONDITIONS 

Appropriate  boundary  conditions  at  the  body  surface  and  the  shock  must 
be  specified  for  the  set  of  governing  equations.  At  the  body  surface  or  wall, 
no-slip  and  no  temperature  jump  conditions  are  used.  Thus,  u^,  =  v^,  =  =  0, 

and  the  wall  temperature  or  heat-transfer  rate  is  specified.  The  conditions 
immediately  behind  the  shock  are  obtained  from  the  Rankine-Hugoniot  relations. 

EDDY-VISCOSITY  MODEL 

For  turbulent  flow  a  two-layer  eddy-vitcosi  model  introduced  by  Cebeci® 
consisting  of  an  inner  law  based  upon  Prandtl 's  mixing-length  concept  and  the 
Klebanoff^'-Clauser^^  expression  for  the  outer  law  is  used.  Further  details 
on  the  eddy-viscosity  model  can  be  found  in  Reference  5. 


THERMODYNAMIC  AND  TRANSPORT  PROPERTIES 

For  a  perfect  gas,  the  thermodynamic  properties  for  specific  heat  and 
enthalpy  can  be  expressed  as 

Cp  =  yR/(y-i) 

h  -  C  T 
P 
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The  viscosity  is  calculated  from  Sutherland's  viscosity  law: 

3/2 

y  =  2.27  E-8  T  '  /(T  +  198.6)  slug/ft-sec 

The  Prandtl  number  is  assumed  constant  everywhere. 

For  air  in  chemical  equilibrium  a  table  look-up  procedure  or  curve-fit 
data  is  used  to  provide  the  thermodynamic  and  transport  properties  as  a  function 
of  the  pressure  and  enthalpy. 

TABLE  LOOK-UP 

A  two-dimensional  table  was  generated  for  the  properties  using  the  method 
developed  by  Miner,  Anderson  and  Lewis. ^  For  a  given  pressure  and  temperature, 
the  enthalpy  and  density  are  determined  using  the  reservoir  calculations  of 
Lordi.^**  The  viscosity  is  obtained  by  curve-fits  from  the  Wilke  semi-emperical 
formula,^  while  the  Prandtl  number  is  obtained  by  interpolation  of  the  Hansen 
data.® 

COHEN  CURVE-FITS 

Curve-fit  data  are  based  on  Cohen's  fit^  of  Hansen's  tables®  for  the 
transport  properties  and  Moeckel  tables'®  for  the  thermodynamic  properties  of 
equilibrium  air. 

Density:  The  enthalpy  dependence  of  the  density  is  given  by  the  curve-fit: 


—  =  1.0 

P 


-  1.0477 


r  /  .  \0.6123T 

r-(^)  ] 


This  fit  is  reasonably  good  for  the  enthalpy  range  0.0152  <  h/h  2.0. 


The  maximum  deviation  in  this  range  is  about  25  percent  at  low  enthalpy,  and 
the  average  deviation  for  all  data  is  about  six  percent. 


The  pressure  dependence  is  given  by  the  following: 

965 


0.0294  I  — 


(%)•' 


This  fit  has  a  deviation  of  less  than  one-half  percent  over  the  range  of 
pressures  10“^  £  p/p-  £  10. 


Viscosity:  The  enthalpy  dependence  for  the  viscosity-density  product  is  given 
by  the  curve-fit: 


Pe^e 


P  /  \  0.3329  -I 


PV 


1.0  -  1.0213 


This  fit  has  better  agreement  than  that  for  the  density.  The  maximum 
deviation  over  the  enthalpy  range  0.0152  ^  —  2.0  is  about  8  percent  and 

the  average  deviation  is  about  three  percent. 


For  the  pressure  dependence  the  following  was  used: 


1 

I 

I 

1 

1 

1 

1 


Pc^c 


=  0.225 


{kr 


This  fit  too  has  a  deviation  of  about  one-half  percent. 


Prandtl  Number:  The  dependence  of  Prandtl  number  on  pressure  is  neglected, 
and  its  variation  with  enthalpy  is  assumed  to  be  of  the  form  shown  in  Figure 
2,  which  is  a  fit  of  Hansen's  data.®  A  pre-condition  for  the  validity  of  an 
effective  Prandtl  number  is  that  the  gas  must  be  in  local  thermodynamic 
equilibrium. 

BADE  CURVE-FITS 

The  density  curve-fits  from  Cohen^  are  not  in  good  agreement  at  low 
enthalpy.  Better  values  of  density  are  obtained  using  curve-fits  suggested  by 
Bade“ 

k  ■  ik) ik) 

where, 

X  =  0.70  +  0.04  logj^Q  31. 9<  h/RT^  _<  480 

X  =  0.94  if  h/RT^  <  31.9 


The  fits  have  a  maximum  deviation  of  seven  percent. 

The  viscosity  and  Prandtl  number  are  then  evaluated  as  in  the  Cohen 
curve-fits . 

REFERENCE  QUANTITIES  IN  CURVE-FITS 

The  following  reference  quantities  were  used  in  the  calculations: 


h,,  -  2.119  X  10®  ft^/sec^  ■ 

L 

Pj,  =  1.0  atm 

p  =  3.584  X  10  ^  slugs/ft 
«  2.498  X  10  ®  slugs/ft 
pg  -  0.01  Ib/ft® 
hg  -  1080  Btu/lb  -  31.9  RT, 


8465  Btu/lbm 


sec 


1.0  atm 


METHOD  OF  SOLUTION 


Davis^’  presented  an  Implicit  finite-difference  method  to  solve  the 
viscous  shock-layer  equations  for  axially  symmetric  flows.  Murray  and  Lewis'* 
extended  the  method  of  solution  to  three-dimensional,  high  angle  of  attack 
conditions.  The  present  method  of  solution  is  identical  to  that  of  Murray 
and  Lewis.  Therefore,  only  an  overview  of  the  solution  procedure  is  presented 
here. 


The  equations  are  written  in  the  standard  form 


.  .  3W 

*0  ■"*!  r 

9  ri  9  n 


+  AW  +  A  +  A,  —  +  A  —  =0 
2  ^  9s  5  94> 


The  derivatives  are  evaluated  by  the  finite-difference  expressions  used  by 
Frieders  and  Lewis'®  and  substituted  into  the  parabolic  equation  giving  the 
standard  finite-difference  form.  The  difference  equation  can  be  solved  by 
the  method  developed  by  Rlchtmyer. 

The  continuity  and  normal  momentum  equations  are  solved  by  a  similar 
method,  but  they  are  coupled  together.  Finally,  the  shock-standoff  distance 
is  evaluated  by  integrating  the  continuity  equation  as  discussed  by  Murray 
and  Lewis. '* 


The  solution  begins  on  the  spherically-blunted  nose  by  obtaining  an  axi- 
symmetric  solution  in  the  wind-fixed  coordinate  system.  At  a  specified  location, 
the  axisymmetrlc  solution  is  rotated  into  the  body-fixed  coordinates  and  used 
as  the  initial  profile  for  the  three-dimensional  solution.  The  three- 
dimensional  solution  begins  on  the  windward  plane  and  marches  around  the  body 
obtaining  a  converged  solution  at  each  <{>  step.  After  completing  a  sweep  in 
<(),  the  procedure  then  steps  downstream  in  s  and  begins  the  next  (J)  sweep. 

At  each  point  the  equations  are  solved  in  the  following  order:  (i)  (J)-momentum, 
(ii)  energy,  (iii)  s-momentum,  (iv)  integration  of  continuity  for  ngjj  and  (v) 
the  coupled  continuity  and  normal  momentum  equations. 


RESULTS  AND  DISCUSSION 

The  freestream  conditions  used  for  the  three  altitudes  are  tabulated  in 
Table  1,  while  the  sample  cases  presented  are  shown  in  Table  2.  The  flow  was 
assumed  to  be  turbulent  at  AO, 000  feet  and  laminar  at  80,000  and  160,000  feet 
altitudes.  The  dimensions  of  the  seven-degree  half-angle  spherically  blunted 
cone  are  shown  in  Figure  1.  The  angles  of  attack  considered  were  0,  2,  10  and 
20  degrees.  The  three  wall  conditions  were  a  cold  wall  (5A0R) ,  a  cool  wall 
(3600R)  and  an  adiabatic  wall. 

TEST  CASE  RESULTS 

A  wide  range  of  conditions  was  attempted,  but  since  the  variations  were 
similar,  only  a  few  sample  results  are  presented. 

There  is  a  significant  change  in  the  shock  standoff  distance  for  perfect 
gas  and  air  in  chemical  equilibrium  (Figure  3) .  The  shock-layer  thickness  is 


nearly  the  same  for  all  options  of  the  equilibrium  properties. 

Case  A  (a  =  0  degrees): 

The  streamwise  variation  of  wall  pressure  with  altitude  shows  very  little 
differences  (Figure  4)  at  AO, 000  and  80,000  feet  altitudes,  but  lower  values 
are  predicted  for  160,000  feet.  These  values  are  nearly  identical  for  all  the 
gas  models  except  within  the  nose  region. 

The  streamwise  variation  of  wall  heat-transfer  rate  is  shown  in  Figure  5 
for  all  the  cases.  In  general,  the  perfect  gas  and  table  look-up  values  are 
the  lowest  and  highest  respectively,  with  the  Cohen  and  Bade  values  lying  in 
between. 

The  variation  of  streamwise  skin-friction  coefficient  shows  trends  similar 
to  those  for  the  wall  heat-transfer  rate.  For  some  cases,  the  Cohen  and  Bade 
results  do  not  lie  within  the  band  of  perfect  gas  and  table  look-up  values. 

The  force  and  moment  data  (Table  3)  show  a  nine  percent  difference  in  the 
axial  force  coefficients  from  the  perfect  gas  and  table  look-up  results,  while 
these  differences  are  quite  small  for  the  Cohen  and  Bade  results. 

Case  B  (cx  =  2  degrees) 

The  variation  of  the  wall  pressure  for  the  windward  plane  is  shown  in 
Figure  7.  The  pressures  from  the  table  look-up  and  Bade  options  are  nearly 
identical.  On  the  spherical  nose,  the  differences  between  perfect  gas  and 
table  look-up  are  substantial,  and  after  tending  together  for  some  distance, 
they  again  tend  to  diverge. 

The  wall  heat-transfer  rate  for  the  windward  plane  from  the  perfect  gas 
prediction  is  about  ten  percent  less  than  the  table  look-up  results  (Figure  8) 
and  diverges  downstream.  The  difference  between  the  table  look-up  and  Bade 
data  is  fairly  small  (within  two  percent) . 

The  streamwise  skin-friction  coefficient  shows  results  similar  to  the 
wall  heat-transfer  rate  (Figure  9)  data. 

The  transverse  variation  of  wall  pressure  (Figure  10),  wall  heat-transfer 
rate  (Figure  11),  streamwise  (Figure  12)  and  transverse  (Figure  13)  skin- 
friction  coefficients  show  similar  trends. 

The  variations  in  the  force  and  moment  data  for  one  case  are  shown  in 
Table  A.  The  inviscld  axial  and  normal  force  coefficients  for  perfect  gas 
conditions  are  four  and  twelve  percent  less  than  the  table  look-up  results, 
while  the  viscous  axial  and  normal  force  coefficients  are  nine  and  sixteen 
percent  less.  The  variation  of  the  center  of  pressure  is  about  two  percent. 

The  differences  between  the  table  look-up  and  Bade  results  are  much  smaller 
and  quite  acceptable. 

Case  C  (a  •  10  degrees) 

The  variation  of  wall  pressure  for  the  windward  plane  does  not  show  the 
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divergence  noted  earlier,  but  substantial  differences  are  observed  in  the  nose 
region  (Figure  14). 

The  streamwise  skin-friction  coefficient  (Figure  15)  shows  good  agree¬ 
ment  for  this  case. 

The  transverse  variation  of  wall  pressure  (Figure  16)  shows  large  dif¬ 
ferences  on  the  leeward  side.  Due  to  crossflow  separation,  data  are  not 
available  up  to  the  leeward  plane.  The  differences  in  streamwise  skin-friction 
coefficient  are  nearly  constant  in  the  transverse  direction  at  a  downstream 
section  (Figure  17).  The  transverse  variation  of  the  transverse  skin-friction 
coefficient  is  shown  in  Figure  18. 

The  force  and  moment  data  for  this  case  are  shown  in  Table  5.  The 
differences  between  perfect  gas  and  the  table  look-up  results  are  much  smaller. 
The  axial  force  coefficient  for  perfect  gas  conditions  is  seven  percent  less 
than  the  table  look-up  results,  while  the  center  of  pressure  locations  are 
within  about  one  percent. 

Case  D  (a  =  20  degrees) 

The  streamwise  variations  for  all  the  quantities  are  similar  to  those  for 
Case  C.  At  the  same  section  there  is  earlier  crossflow  separation  due  to  the 
higher  angle  of  attack. 

The  transverse  variations  of  wall  pressure  (Figure  19),  wall  heat-transfer 
rate  (Figure  20),  streamwise  (Figure  21)  and  transverse  (Figure  22)  skin- 
friction  coefficients  are  also  similar. 

The  axial  force  coefficient  for  perfect  gas  is  six  percent  less  than  the 
table  look-up  results  (Table  6),  while  the  center  of  pressure  location  is 
nearly  the  same. 

Computing  Times  Required 

Computing  times  given  in  the  tables  are  with  the  table  look-up  time  as  a 
reference  of  100.  Due  to  the  absence  of  some  input/output  for  the  perfect  gas 
model,  the  time  for  this  condition  is  slightly  lower  than  the  Cohen  and  Bade 
options.  If  one  omits  the  extra  printout,  the  curve-fit  of  equilibrium  proper¬ 
ties  and  perfect  gas  computing  times  are  very  comparable  and  result  in  a  sixty 
percent  reduction  in  the  computing  time  required  for  the  table  look-up  model. 

For  low  values  of  wall  temperature,  since  all  the  equilibrium  models  es¬ 
sentially  use  Sutherland's  viscosity  and  the  density  from  the  table  look-up 
and  Bade  options  is  in  good  agreement,  the  differences  in  the  wall  heat- 
transfer  rate  are  probably  due  to  the  differences  in  the  Prandtl  number. 


CONCLUSIONS 

In  general,  wall  heat-transfer  rate  differences  between  the  perfect  gas 
and  table  look-up  results  are  about  ten  percent,  while  the  axial  force  coef¬ 
ficients  differ  by  five  to  fifteen  percent.  Large  pressure  differences  in  the 
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nose  region  (which  constitutes  a  significant  percentage  of  the  frontal  area) 
probably  cause  the  large  differences  in  axial  force. 

Curve-fits  can  be  effectively  used  to  model  the  thermodynamic  and  trans¬ 
port  properties  of  air  in  chemical  equilibrium  producing  reasonably  good  force 
and  moment  coefficient  predictions  and  with  computing  times  that  are  com¬ 
parable  to  the  perfect  gas  model.  The  differences  in  wall  heat-transfer  rate 
for  equilibrium  flow  are  probably  caused  by  the  differences  in  Prandtl  number. 

Since  there  are  considerable  differences  between  the  perfect  gas  and 
equilibrium  air  results,  the  viscous  shock-layer  method  just  developed  should 
be  used  to  predict  viscous  flows  in  chemical  equilibrium,  such  as  for  reentry 
conditions,  with  good  accuracy  and  within  reasonable  computing  times. 


NOMENCLATURE 

skin-fricticn  coefficient  in  the  streamwise  direction 
skin-friction  coefficient  in  the  transverse  direction 


constant  pressure  specific  heat 

2 

static  enthalpy,  h*/Uoo 

2  2  2 

total  enthalpy,  h  +  (u  +  v  +  w  )/2 

- 2 

eddy  thermal  conductivity,  -{pv)'h'/(e  9h/8n) 


Mach  number 

shock-standoff  distance,  n  , /R* 

sh 

pressure,  p  /p  U 


Prandtl  number 


turbulent  Prandtl  number,  C  N^/k^ 

2 

convective  heating  rate,  Btu/ft  -sec 

A  A 

distance  from  and  normal  to  the  body  axis,  r  /R 
body  nose  radius 


general  surface-normal  coordinate  system  streamwise,  normal  and 
crossflow  directions 

T  temperature 

it  2  * 

reference  temperature,  U^/C^ 
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u,v,w  streamwise,  normal  and  crossflow  velocity  components 
by  the  freest  ream  velocity,  U* 

W  dependent  variable 

a  angle  of  attack 

y  ratio  of  specific  heat 

e  Reynolds  number  parameter,  r  =  u  /n  U  R 

^ref  *'oo  00 

6  body  angle  in  the  streamwise  direction 

•k  k 

y  viscosity,  y 

yj.g£  reference  viscosity,  y*{T 

viscosity,  - (pv)"u'*/ (eSu/8n) 
p  density,  p*/p^ 

Subscripts 

Bade 
Cohen 
Inviscid 
Streamwise 

conditions  behind  the  bow  shock  wave 
turbulent  quantity 
wall  condition 

dimensional  freestream  conditions 
transverse 
Superscripts 

* 


B 

C 

INV 

s 

sh 

t 

w 

00 

4* 


dimensional  quantity 
fluctuating  value 
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TABLE  1.  Freestream 

Conditions 

Altitude,  h(feet) 

40,000 

80,000 

M 

00 

20.6 

20.4 

Velocity,  U^j^Cft/sec) 

20,000 

20,000 

2 

Pressure,  p^(lb/ft  ) 

393.085 

58.505 

Temperature,  T^(R) 

389.97 

397.69 

3 

Density,  (slug/ft  ) 

5.873E-4 

8.571E-5 

160,000 

18,4 

20,000 

1.942 

487,17 

2.322E-6 


Reynolds  Number  (per  foot) 


3.955E7 


5.679E6 


1.305E5 


TABLE  2.  Test  Cases 


cc 

(degrees) 

Alt 

(feet) 

T 

w 

(R) 

T  /T 
w  “ 

Flow 

0 

80,000 

540 

1.36 

Laminar 

2 

40,000 

3600 

9.16 

Turbulent 

10 

80,000 

Adiabatic 

- 

Laminar 

20 

80,000 

540 

1.36 

Laminar 
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TABLE  3.  Force  and  Moment  Data  (Case  A) 


a  =  3  deg,  Alt  =  80,000  feet,  =  540R,  S/R  =  30 


Actual  Values 


PC 

TLU 

COHEN 

'‘lOV 

0.070216 

0.076465 

0.076465 

% 

0.06831 

0.07538 

0.07509 

0.00251 

0.00245 

0.00209 

Ca 

TOT 

0.070824 

0.077835 

0.077176 

TIME* 

(Sec) 

72 

233 

Percent  Variation  Based 

102 

on  TLU 

PC 

COHEN 

BADE 

8.2 

0.0 

0.0 

% 

9.4 

0.4 

-0.3 

-2.4 

14.7 

11.0 

^OT 

9.0 

0.8 

0.1 

time’'" 

33 

44 

43 

*IBM  370/3032  with  FORTHX  COMPILER  (OPT2) 

^Based  on  TLU=100 
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BADE 

0.076465 

0.07557 

0.00218 

0.077757 

100 


■ 


I 


TABLE  4.  Force  and  Moment  Data  (Case  B) 


a  =  2 

deg,  Alt  =  40,000 

feet,  T  =  3600°R  S/R 

w 

=  32 

Percent 

Actual  Values 

Variation 

Based  on 

TLU 

PG 

TLU 

BADE 

PG 

BADE 

C* 

0.06331120 

0.0654952 

0.0654952 

3.3 

0.0 

^INV 

"^INV 

0.0345556 

0.0392156 

0.0392156 

11.9 

0.0 

"a 

0.071163 

0.078564 

0.078414 

9.4 

0.2 

s 

0.032386 

0.038607 

0.038653 

16.1 

-0.1 

s 

-0.021207 

-0.025800 

-0.025856 

17.8 

-0.2 

Z  /L 
cp 

0.654814 

0.668288 

0.668924 

2.0 

-0.1 

'*SF 

0.00595 

0.00744 

0.00710 

20.0 

4.6 

0.06521 

0.07112 

0.07132 

8.3 

-0.3 

TIME 

1778* 

5203* 

2384* 

34^ 

46' 

* 

Seconds,  IBM  370/3032  with  FORTHX  COMPILER  (0PT2) 
^Based  on  TLU=100 
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TABLE  5.  Force  and  Moment  Data  (Case  C) 


a  =  10  deg,  Alt  =  80,000  ft,  =  Adiabatic,  S/R  =  32 


Actual  Values 

Percent  Variation 
Based  on  TLU 

PG 

TLU 

PG 

INV 

0.094359 

0.097347 

3.1 

INV 

0.296715 

0.302710 

2.0 

^A 

0.099410 

0.106954 

7.1 

S 

0.302246 

0.299984 

-0.8 

S 

-0.198296 

-0.194445 

-2.0 

Z  /L 
cp 

0.656076 

0.648185 

-1.2 

0.00433 

0.00482 

10.2 

C 


0.09508 


0.10214 


6.9 


TABLE  6.  Force  and  Moment  Data  (Case  D) 


20  deg,  Alt  =  80,000  ft,  =  5A0R,  S/R  =  32 


Actual  Values 

PC, 

TLU 

Percent  V'ariation 
Based  on  TLU 

PG 

0.151591 

0.155653 

2.6 

0.716164 

0.727595 

1.6 

0.159610 

0.170116 

6.2 

0.719390 

0.732410 

1.8 

-0.450558 

-0.459907 

2.0 

0.626305 

0.627937 

0.3 

0.00613 

0.00644 

4.8 

0.15348 


0.16368 


6.2 


Flgur 


Figure  2.  Prandtl  Number  Variation  with  Enthalpy  Ratio  (Cohen  and  Bade  Options) 


Streamwise  Variation  of  Shock  Standoff  Distance 


Streamwlse  Variation  of  Wall  Heat-Transfer  Rate  (Case 


Streamwise  Variation  of  C,  (Case  A.  a  =  0  6m.  Alt  =  An  nnn 


Streamwise  Variation  of  Wall  Pressure  for  <t>  =  0  Plane 


Streamwise  Variation  of  Cf  for  $  =  0  Plane  (Case 


o 


1+3’ I  0+3‘l 

dNId/Md 
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Figure  10.  Transverse  Variation  of  Wall  Pressure  at  S/R  = 
(Case  B,  a  =  2  deg,  Alt  =  40,000  ft,  T  =  3600R) 


Transverse  Variation  of  C.  at  S/R  =  20 

(Case  B,  a  =  2  deg,  Alt  40,000  ft,  T  -  3600R] 


Transverse  Variation  of  Cf .  at  S/R  =  20 

(Case  U,  u  =  2  deg,  Alt  -  AO, 000  ft,  T  ■»  3600R) 


Streamwise  Variation  of  Cf  for  d)  =  0  Plane 


Figure  16.  Transverse  Variation  of  Wall  Pressure  at  S/R  =  20 

(Case  C,  a  =  10  deg,  Alt  =  80,000  ft,  T  =  Adiabatic) 


Figure  17.  Transverse  Variation  of  Cf^  at  S/R  = 
(CuHi-  C;,  a  -  iO  Alt  -  MO.ODU  It. 


Transverse  Variat 


B7.5G 


